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AEROELASTIC STUDY OF A SIMPLE MODEL 



ABSTRACT 

This thesis concerns the study of a two dimensional model subjected to 
gust loads. The model consists of a wing attached to a fuselage of infinite 
moment of inertia by means of a tension-compression spring and a torsion 
spring. 

The equations of motion were developed using two-dimensional thin airfoil 
theory for incompressible flow and they are solved on an electronic analogue 
computer. Although the equations are set up to handle a gust of any arbitrary 
gradient, the analysis was made for a sharp-edged gust and comparison then 
made with gusts of finite gradient. 

It was found that stresses in the wing were reduced by locating the center 
of gravity near the elastic center and by reducing the radius of gyration of the 
wing. There were two oscillatory modes, one primarily ben ding and the other 
primarily torsion. The frequencies of oscillation could be lowered by locating the 
center of gravity near the elastic center, by increasing the radius of gyration 
or by increasing the wing to fuselage mass ratio. 

The mode primarily characterized by torsion had a relatively high frequency 
at low flight velocity. The two frequencies approached each other as flight 
velocity increased and damping was improved in each mode until flutter was 
imminent, when the lower frequency (bending) mode very suddenly became un- 
stable. It was found that neglecting the unsteady effects yielded a highly con- 
servative estimate of flutter speed in this particular case. 

Increase in altitude has little effect but to lessen the damping in each mode. 
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This investigation was performed by Lt. Raymond L. Miller USN at the 



University of Michigan in 1955-6. It was done under the supervision of 



Associate Professor Gabriel Isakson. 



AEROELASTIC STUDY OF A SIMPLE MODEL 



INTRODUCTION 

The purpose of the present thesis was the investigation of the response to 
gust loads of a two-dimensional model with three degrees of freedom; particu- 
larly, the effect upon such response of changes in environmental and model 
parameters. 

The present thesis was submitted at the University of Michigan in partial 
fulfillment of the requirements for the degree of Master of Science in Aero- 
nautical Engineering. The subject was suggested by Associate Professor 
Gabriel Isakson of the Aeronautical Engineering Department and performed 
under his supervision. The approach to the problem was patterned after an 
earlier report by Professor Isakson (Ref. 1), but, with solution of the equations 
performed by means of electronic analogue computers, rather than analytically. 

The problem to be investigated is that of aeroelastic effects on a two- 
dimensional wing-fuselage combination in incompressible flow. The wing has 
freedom in twist and bending; the fuselage is free to plunge vertically, but is 



restrained in pitch. 



. 
















2 



THEORY 



1 . The Problem 

The two-dimensional model consists of a fuselage and thin wing of constant 
cross section connected by a tension-compression spring and a torsion spring, 
representing the elastic constraints of the wing in bending and torsion re- 
spectively, as shown in Figure 1. Z w and Zf are vertical displacements of 
the wing and fuselage respectively, relative to fixed axes, and 0 is the 
rotation of the wing chord about the elastic center. All are positive in the 
directions shown and are measured from the initial steady state position. 
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The angles « 0 , 0 and the changes in angle of attack due to gusts and due 
to vertical velocity of the wing Z are assumed to be small and the cosines of 
these angles are taken to be unity in evaluating aerodynamic lift and moment. 
The gust acting upon the model is taken to be normal to the velocity vector, U, 
of the undisturbed free stream. 

2. The Equations of Motion 

The equations of motion of the system are: 



m 



w 



Z w + R0 + k(Z w - Z f ) = -L 



RZ W + I e + K 9 - M 



m f 7^ - k(Z w - Zf) = 0 



( 1 ) 

( 2 ) 

(3) 



where: 



k ; spring constant of tension-compression spring. 

K = spring constant of torsion spring. 

I = mass moment of inertia of airfoil section about the elastic center. 

- ^cg "** ®wO cg C> 

I C g — mass moment of inertia of airfoil section about the center of gravity. 

L r instantaneous value of aerodynamic lift per unit span, as a| departure 
from the initial steady state value, positive up. 

M = instantaneous value of aerodynamic moment per unit span, as a departure 
from the initial steady state value, positive in the direction of 0. 
m£ r mass of fuselage per unit span 
m^ s mass of wing per unit span 
R = static mass moment about the elastic center. 



= “wflcg 0 ) 
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It is convenient to change from a time variable, t, to a variable, s, 
where: 

s = U t = distance travelled in half-chords. 

_ 2 

U = U^; = velocity of flight in half-chords per unit of time. 
Under this transformation: 



d 
dT 

d2 

dP = 



_ d 

*-ar 



= u 



^ 2 d 2 



Denoting differentiation with respect to s by prime marks: 
2 

m w U Zw + RU 2 6” + k(Z w - Z f ) = -L<s) 
RU 2 Z^ + IU 2 6 " + K6 = M(s) 
m f U 2 Z’^ - k(Z w - Z f ) - 0 
These equations are now put in non-dimensional form: 



(4) 

(5) 

( 6 ) 



7 n 



, /n M Zw Zf 

c + *cg ^ + K B*~Cn " ~C~) - 



w 



L * 



m w U 2 C 



(7) 



l cg -~- 1- r 2 0” - K 0 ^ 

c e C T m u l; 



( 8 ) 



m f Zf” 

m w C 



- 



J w 



B' C 



-) = o 



(9) 



where: 



r “ ~7T y J m~ ~ ra< ^ ius °f gyration as a fraction of the chord 
k kC 2 

Kg - — — ^ = bending stiffness parameter 



„ K K 

= — — {j 2 q2 - 4 ^ ~~ t; 2 — - torsion stiffness parameter 






I 
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3. Aerodynamic Force and Moment 

Four types of force and moment act on the airfoil as a departure from 
initial steady state values: 

1) Gust effect 

2) Wake effect 

3) ’’Quasi-steady'' forces - forces which would be produced by the 
instantaneous velocity of the airfoil in the absence of wake effect. 

4) Apparent mass effect 

The last three forces listed are results of the disturbed motion of the 
airfoil. The quasi-steady forces and the wake effect can be accounted for 
simultaneously; thus, lift and moment each consist of three components: 



L = 


L g + 


L m 


+ Lma 


( 10) 


M = 


M + 

S 


M 

m 


+ 

ma 


(11) 



Lg = lift due to gust effect. 

L m = lift due to quasi-steady forces and wake effect resulting from 
airfoil motion. 

L ma = due to apparent mass effect. 

M z moment due to gust effect. 

S 

M m = moment due to quasi-steady forces and wake effect resulting from 
airfoil motion. 

M ma = moment due to apparent mass effect. 

Lift and moment are computed by thin airfoil theory. 

3.1. Gust effect 

For sharp-edged gusts, the change in lift coefficient is a function of time 



(or distance in half-chords): 



1 



6 



C 



L 



dC L 

d* 



A cx^s) - 2 /T^ — YMs) 



where w is the gust velocity, positive up, and ^(s) is the Kussner function 
for build-up of lift due to change of angle of attack in a sharp-edged gust, the 
value of this function approaching unity with increase is s. 

A gust of arbitrary gradient can be considered to be made up of a series of 
infinitesimal sharp-edged gusts as shown in Fig. 2. The lift coefficient can be 
computed for each infinitesimal gust by the formula above and all the values 
summed by integration. Thus the change of lift coefficient becomes a function 




of (s - cr), where cr is the value of s at the time of origin of the particular 
infinitesimal gust (Aw): 

C L = ^ w W[s - cr) 



Summing the effects of all the infinitesimal gusts: 



2 rr 



9 

A w (<r) Y/(b - or) 



} 



U 



w(o) Y' (s) + 



C r 



2 TT 
U 
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w(o)Yms) - O) Act 



CTr 



or, in the limit as Acr -*■(): 
2T7 



p- |w(o) y/ (s) 4 Jw-(cr) y/ (s - cr > c/<r| 



( 12 ) 



This is the Duhamel, or superposition, integral. An equivalent form is found 
by integration by parts: 



4f— | w (s)/'(o) ■+ j S w (<r) y'(a -<r) dcr^j 



(13) 



It will be shown later that Y'(O) = 0, thus: 



C, r 277 



w (cr) 



— - V y’(s-cr) 

u r 



dcr 



L g = -f cu2 c l ■ -r c3u2c l 



_ j W 

L g (s) = — U 2 C 3 J — (cr) £r)dcr 

O 



(14) 



In two-dimensional thin airfoil theory for incompressible flow, the aero- 
dynamic center is at the 25% chord point: 

*ac = ^25 

M g (s) = -(l ac C)L g (s) (15) 



3. 2 Quasi-Steady Forces and Wake Effect 

The sum of the quasi-steady forces and wake effect is a function not only of 
the instantaneous disturbed motion of the airfoil, but also of the time history of 
that motion. These forces build up with time in much the same fashion as 
change in forces due to a gust and are handled in similar fashion. It can be 
shown that the terms which are functions of the time history of the disturbed 
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motion are all related to the quantity: 

<z w -+ U0 +1 75 C9) 

(R e£. 2, page 282). This quantity can be seen to be the instantaneous vertical 
velocity of the three-quarter chord point and is equal to the angle of attack at 
that point times the flight velocity, U. Thus, by resolving velocities at the three- 
quarter chord point, the quasi-steady forces and wake effect can be accounted 
for simultaneously The lift due to these forces can be found by a process 
analogous to that due to the gust effect: 

/■ ^ 

L m = U 2 C 3 j «^ 5 (<r)dcr (16) 

O 

where: 

_ I 75 C 1 

c <75< t > = 9<‘> + — [j— 9 (t) + Z w< 1 > 

Z' 

<* 75 (s) = 0(s) + 21 ?5 0 '(s) + 2 C W ~ (17) 

o( is the angle of attack at the three-quarter chord point, measured as a 
7 o 

departure from the initial steady state value, and ^s) is the Wagner function 
for build up of lift coefficient following a sudden change of angle of attack. Like 
the Kussner function, the Wagner function is defined to approach unity with 
increasing s. 

The above equation is a form of the Duhamel integral as in equation (12), 
but with ( 0 ) - 0 . 

The lift due to the quasi-steady forces and the wake effect acts at the 25% 

C 2 

chord point, except for that part of the quasi-steady lift, equal to 77/° (— ^ — ) (U 8 ), 
which is due to the rotational velocity of the airfoil and acts at the 50% chord 
point (Ref,. 3, page 104). This, then, must be subtracted from the moment: 



c 
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TTP ^~) 2 -^-0’(s) z U 2 C 4 Q'( s) 



M m (s) = 



-‘ac L m - ^U 2 C 4 0'( S ) 



( 18 ) 



3.3. Apparent Maas 

The apparent mass lift due to the disturbed motion of the airfoil may be 

thought of as the force associated with the acceleration of a cylinder of air 

C 2 

with diameter equal to the wing chord. It is equal to the product of TTP { — ) 

and the instantaneous acceleration of the 50% chord point. In addition a pure 

TTP C 4 

couple is acting on the airfoil equal to the product of — g — (-j— ) and the 
instantaneous rotational acceleration of the airfoil (Ref. 2, page 263). 

The distrubed velocity of the 50% chord point relative to the free stream is: 

z w + U0ti 5o ce 

and the acceleration of the point is: 



Z w + U0+l 5O C0 



Thus; 



ma 



77P o •• 

~r~ c 2 (Z w •+■ us+i 5 O c0) 



TT^ D 2 C 3 ( 



Z w 1 

IT + TF e ’ * ‘so*"> 



(19) 



M = - Ica CL 
ma 50 ma 



8 16 



- - 1 CL 

ac ma 



- Tr- 'so* F» e "]< 20 > 



summing equations 14, 16, 19 and 15, 18, 20 for total lift and moment: 



' 
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L(s) r ^-U 2 C 3 j j^-(cr) ^’(s -cr )dc^ -tr )0(^(cr)d<r 



+ (4 1 + 4- 0' + 1 6O 0 m ) 



( 21 ) 



M <s) = -1^ C L(s) " 



TTP -■ 



2r>4 



16 



U^C 



“cr 



+ (1 50 + ~ «) 0 " ( (22) 



2 -2 2 
Dividing equation (21) by m D C and equation (22) by c^U C 



w 

$ 



L ( s ) 1 I ( w 

m w U 2 C = ~^T j)lT (<r) y {B " <r > dcr + j^( s - cr )^ 5 (0') d£r 



+ &L 



c * ~ 0 + *50 Q M) 



(23) 



M(s) 

m U 2 C 2 
w 



1 L(s) 
A ac 

m U 2 C 
w 



\JJ / c 






- - l^Vi* ~^) d ^ + /^(s -^^5 ( cr )dcr 



2* 



- — /] 

/y / J 50 c 



75 



*’ * < 1 50 + "^~>a”}(24) 



where: 



M - - ,r\r^ - relative density 

= ratio of mass of wing to mass of a cylinder of air with 
diameter equal to wing chord. 



Now these expressions for lift and moment can be substituted into the equations 



of motion, (7) and (8): 



11 



+ < leg * 7T ) 9" + Tp Q' +k b( 



2vv; 

C 



_ Zf 



fj V'** 3 ~ °~) d<r + j <P(-S-<r) : *ir( a ')J<rj - O (<?5) 



tr* + 



~^)e" + n cz .!g)Z + ig e’+ « r e 

Jj* f f ft ^ X* ^S'CT/c/o- -f f 0 Cs -cr^ ^ 7 y ^ ) J<r^j ~ O (26) 



4. Wagner and Kussner functions 

The exact forms of the Wagner and Kussner functions for incompressible 
flow involve Bessel functions (Ref. 2, page 285, 287): 



where: 




F(k) + , G(k) = C (k) = 



(2) 

H n = Jn “ iY n 



^^(k) , H 0 ^'(k) 



(2)/ 



(k) 



F G (k) + i G G (k) = C(k) [j 0 (k) - IJ^k)] ■* ,J X 

J n = Bessel function of the first kind of order n 

Y = Bessel function of the second kind of order n 
n 

Since these functions are too complex for analytical use, it is necessary to use 
the following approximations: 

0(s) - 1 - 0.165 e '°- 0455s - 0. 335 e~°* 35 (27) 

0(3) = 1 - 0.500 e*°* 130s - 0„ 500 e _s (28) 



These equations show that a sudden change in angle of attack due to an instan- 
taneous disturbed motion of the airfoil immediately causes a change in lift 
coefficient then increases asymptotically to the ultimate steady state value; 
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whereas, a sharp-edged gust causes no initial change in lift coefficient, but then 
causes the lift coefficient to increase asymptotically to its ultimate steady state 
value. 

Using the approximate equations, (27) and (28), for the Wagner and Kussner 
functions, the equations of motion can be solved, either analytically or by means 
of the analogue computer. 

i 

5. Dynamic Overstress Factors 

The dynamic overstress factor is defined as the ratio of the maximum stress 
the wing of the elastic model is subjected to by a given gust to the maximum stress 
the corresponding rigid model (K-j. = Kg = o ) would be subjected to by the same 
gust. 

Dynamic overstress factor in bending 

. * elastic = <q> „ 

K ( ) rigid (Z't ) rigid 

C C 

Dynamic overstress factor in torsion 

(K Q ) elastic 
( K 9 ) rigid 

From equations (8) and (24): 

(K^B) rigid . m $ c z - 1 

w & 

lgp Zw 
- s c 



I 
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thus; 



<K t 0 ) rigid = 




O 




(29) 



PROCEDURE 



Representative values of the various parameters were selected to serve 
as a standard or basis of comparison in observing gust and flutter effects. The 
standard values of these parameters were selected arbitrarily, but were intended 
to represent a transport or trainer type of aircraft. Standard values and range 
considered are shown in Table I. 



coupled frequencies of approximately 200 cycles per minute in bending and 475 
cycles per minute in torsion: 



The effects of individual model parameters were investigated while holding 
all other parameters fixed. Similarly, the effect of altitude and flight velocity 







eycles/min. 



(jj t = -^5 4 - 470 cycles/min. 



were investigated. Dynamic overstress factors were determined by measuring 
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the maximum acceleration of the fuselage and twist of the wing under a step 
input representing a sharp-edge gust. The period, damping time and mode 
shape of each of the oscillatory modes were determined by adjusting initial 
conditions to separate the modes. 



ELECTRONIC CIRCUITS 



1 . General 



The computer primarily used for solution of the problem was a Reeves 
Electronic Analogue Computer (RE AC Mod. C 101). Traces were made on a 
Brush Oscillograph recorder (Model BL-928) using a Brush Dual Channel D-C 
Amplifier (Model BL-202). 

The complete circuit diagrams are shown in Figs. 5 and 6. 

2. Scaling 

It was found convenient to speed up the solution by scaling. Under the 
transformation, = . 2s, the equations become: 



Z\ 

C 



JL_ = 25 



Zf 



m ^ 



) 



(30) 



*75= 6 + 



. 4 2 &' + o 4r—^~ 

75 c 



(31) 



Z'w 



= -( leg/'- 1 - Iso ) g » - - -aiKiM-V - #*-> 

/< + i y 1 y + 1 



;*/[- 1 1 '*' 



(s - cr)d cr + 



[As- 



(<r) dcr 
/ o 



} 



(32) 



lgo‘3 r)S" = - <l cg /<* 



g/ 



Zft 

— - 4 . Oi 

c 75 



-2SK^Jjd - 251^ j (-w-(cr)y/’ (s -<r)d<r + ( ^(Sj - a )0<^ 5 (<r ) d<r 



(33) 
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$(s,) = 1-0. 165 e - 0. 335 e 



0.2275 B, n 335 e -1.5 s, 
-5 s, „ cnn „-0. 65 s, 



(s, ) - 1 - 0.500 e'° s '- 0.500 e 

A , leg + ^60 //J Z ’^ 1 

(K„ 0) rigid = - ( Lb ' - 



25 



> C 



ac 



3. Symbols 









ft 



~h |^(s - cr )c <j 5 (cr)dcr j 



(34) 

(35) 
’(s -cr )dcr 

(36) 



summer and 
scale changer 

I 



integrator 




■ 0 - 



gain 




resistor (in 
megohms ) 

1 

R 

v/\/\y^— 



I 

capacitor (in 
microfarads) 



amplifier without 
feedback 





1 



amplifier with 1 
megohm feedback 



4, Convolution Integrals 

The analytic solution of the convolution integrals requires an integration 
for each exponential term; similarly the computer solution requires an integrator 
for each exponential term. Consider the following circuit: 
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The differential equation for this circuit is: 

+ ax - - AN(t) 
dt 

and the solution is: 

x = - Ae" at je at N(t)dt 
_ _ Ae" at J e aT N(r) dr 
r - A J < e a<t " 1 ' ) N(T)dr 

O 

Thus the Kussner function can be developed by the circuit shown in Fig. 3. 




Fig. 3. Kussner Function Circuit 



It is possible, by the use of an external circuit, to develop the convolution 
integral with only one amplifier. In order to conserve integrators, the Wagner 
function was handled in such a manner. 
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The above circuit can be shown to satisfy the differential equation: 



} o CRi G o 



e r - 



Rf 

— 5— e! 

R i 1 



and its solution is: 



> _ f* 

o “ R j l 



R ' _i_ (t_>) e i ' (>)d> 



CR i 



and, by the use of helipots, more than one exponential can be handled with 
single amplifier: 

Rf 

<d> 





where: 



R f ( t — ji-rr (t - T) ' 

e ° = ’ A “T ) e e ‘ <T,dr " B 

The circuit for the Wagner function is as shown in Fig. 



i 



t =—(t -r) 

e C 2 R 2 « 




Fig. 4. Wagner Function Circuit 



(r)dr 



<^JS (o') <J<r 



- 
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RESULTS AND DISCUSSION 



1. Accuracy of data 

The accuracy in measurement of amplitudes was less limited by the computer 
than by the recorder, which was accurate within about 4%. This produced good 
results in determination of dynamic overstresses and smooth curves were obtained 
in all cases. Determination of mode shapes was quite difficult, not primarily 
due to any inherent inaccuracy in the computer, but due to the difficulty encountered 
in obtaining the modes in pure form; by averaging over several cycles, consistent 
values were obtained in most cases. (Mode shapes were measured in bending and 
torsion amplitudes only, phase angle changes were too small to be measured 
except near flutter speed. ) The same difficulty was encountered measuring the 
damping of the oscillatory modes and the averaging procedures were somewhat 
less successful. 

2. Behavior of the Model in the Standard Configuration 

In the standard configuration the elastic center is . 35 chord lengths from 
the leading edge, the center of gravity is . 45 chord lengths from the leading 
edge, the radius of gyration is equal to . 25 chord lengths and the wing to fuselage 
mass ratio is .333; the model has a relative density of 50 at sea level with a 
flight velocity of 380 feet per second. The model experiences dynamic overstress 
factors of 1. 90 in bending and 2. 35 in torsion under the action of a sharp-edged 
gust. There are two oscillatory modes, one primarily involving bending and the 
other primarily involving torsion. The "bending" mode has a frequency of 222 
cycles per minute and damps to one-half amplitude in 3. 7 cycles. The torsion 
mode has a frequency of 5S9 cycles per minute and damps to one half amplitude 
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in 5. 7 cycles. The mode shapes are: 



( 



Z w 




bending “ 



3.69 



chord lengths/ radian 



Z 



w _ 



e 



> torsion = “ °* 108 chord lengths/radian 



where the minus sign indicates that the bending and torsional displacements 
are approximately 180° out of phase. 

3. Effect of the Location of the Elastic Center and Center of Gravity 

The effect of the elastic center and center of gravity locations on the 
dynamic overstress factors under the action of a sharp-edged gust is shown in 
Fig. 7; neither has any appreciable effect on the dynamic overstress in bending; 
however, the dynamic overstress in torsion increases as the elastic center and 
the center of gravity move apart. 

Effects of the elastic center and center of gravity location upon frequencies, 
damping and mode shapes are shown in Figs. 8-10; in interpreting these curves, 
it should be remembered that the elastic center is initially at the . 35 chord point 
and the center of gravity at the . 45 chord point and each parameter is adjusted 
individually, not simultaneously. Fig. 8 shows that the frequency of the torsion 
mode is at a minimum when the elastic center is coincident with the center of 
gravity and increases rapidly as the elastic center moves ahead of the center of 
gravity. The frequency of the bending mode is little affected by the elastic center 
or center of gravity location and is not shown in Fig. 8. 

The damping constants obtained in investigating the effect of varying the 
elastic center location were somewhat erratic, but it is apparent from Fig. 9 
that the damping of the torsional mode decreases as the center of gravity moves 
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aft and as the elastic center moves toward the center of gravity. The damping 
of the bending mode is a minimum when the center of gravity is located near the 
. 35 chord point, although the effect on damping of either center of gravity or 
elastic center location is quite small. 

As seen in Fig. 10, the bending mode is free from any influence of torsion 
when the center of gravity is located just slightly ahead of the elastic center, 
and the torsion mode is free from any influence of bending when the elastic 
center is just slightly ahead of the center of gravity. The slight coupling that 
occurs when the elastic center and center of gravity are coincident is due to the 
apparent mass aerodynamic effect. 

4. Effect of the Radius of Gyration 

Effect of the radius of gyration on dynamic overstress factors and mode 
characteristics is plotted in Figs. 11-12. The dynamic overstress factor in 
torsion increases as the radius of gyration increases. This is consistent with 
increases in dynamic overstress factor as the elastic center and center of 
gravity move apart. 

Again in this case, the dampin g constants were somewhat erratic, but damping 
of the torsion mode was seen to improve as the radius of gyration decreased; the 
frequency of the torsion mode increases as the radius of gyration decreases. The 
influence of torsion in the bending mode increases as the radius of gyration becomes 
smaller. The frequency and damping of the bending mode and the mode shape of 
the torsion mode are little affected by the radius of gyration and these effects are 
not shown in Fig. 12. 

5. Effect of Wing to Fuselage Mass Ratio 



The effect of changing the fuselage mass, while holding the wing mass fixed 
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is shown in Figs. 13-15. The dynamic overstress factors are increased some- 

m 

what as the mass ratio, w / m^, is increased. The frequency of the bending 

mode increases and the damping decreases as the fuselage becomes lighter. 

The frequency and damping of the torsion mode is not noticeably affected and is 
not shown. Coupling between the two modes is at a minimum at a mass ratio of 
about . 30. 

6. Effect of Flight Velocity 

The effects of flight velocity on behavior of the model are shown in Figs. 16- 
18. The flutter speed for the model was found to be 1142 feet per second (com- 
pressibility effects neglected). As the flutter speed is approached the frequencies 
of the two modes approach each other as would be expected (Ref. 4, page 222) 
and the damping improves in each mode until very near the flutter speed, where 
the damping of the bending mode drops off sharply. This sudden deterioration of 
damping also might be expected; a study of a three dimensional wing of large 
aspect ratio produced frequency and damping characteristics very similar to 
those of Fig. 17 (Ref. 2, pages 552-3). Torsion becomes more predominant in 
each mode with increase in speed. 

The dynamic overstresses increase gradually at first with increase in flight 
velocity, and then rise sharply at the flutter speed; this is especially true of the 
dynamic overstress in torsion. 

7. Effect of Altitude 

Figs. 19-20 show the effect of altitude or relative density on dynamic over- 
stress factors and damping. As the altitude is increased above sea level (or 
density reduced), it was found that frequencies of the modes are not discemably 
changed, nor are the mode shapes significantly affected; thus frequencies and 
mode shapes are not plotted. The damping is decreased in both modes; this 
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decrease in damping does not mean that the flutter speed would be reduced; in 
fast, it is shown in Ref. 2 (page 542) by analysis of an airfoil model with freedom 
in bending and torsion, that increase in the relative density (above a fairly small 
figure) increases the flutter speed. This is logical since damping is improved 
with increased airspeed. Dynamic overstress factors are not greatly affected 
by the value of the relative density, except for values of relative density 
approaching zero. 

8. Neglect of Unsteady Forces 

Substitution of quasi-steady forces for the unsteady forces is accomplished 
by replacing the Wagner function with unity; this gives a slightly conservative 

i 

estimate of dynamic overstres6es and a very conservative estimate of flutter 
speed - 498 ft/sec. as compared to 1142 ft/sec. It is of interest, however, 
that flutter occurs in the torsion mode, not in the bending mode as with the un- 
steady forces included. An analysis of a jet transport wing (Ref. 2, pages 567-8) 
produced similar results, except that flutter occurs in the torsion mode for both 
unsteady and quasi-steady theory. 

On the basis of a constant derivative analysis, in order for energy to be 
extracted from the air the following criterion must be met. 

WC < 1. 56 
U 

(Ref. 5, page 64). Flutter of this model under quasi-steady analysis occurs at 
a frequency of 58. 7 radians per second; i. e. ; 



U 

9. Gust Gradient 



A gust with a finite gradient would cause lower stresses on the wing than 



I 

I 

I 
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a sharp-edged gust. Fig. 21 shows the effect on the stresses as function of 
time required for the gust to build up to its maximum value. 
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CONCLUSIONS 

It was found that, for this two-dimensional model, dynamic overstresses 
caused by a sharp-edged gust could be reduced by locating the center of gravity 
close to the elastic center or by reducing the radius of gyration of the wing. The 
frequencies of the oscillatory modes were increased by locating the center of 
gravity near the elastic center, by increasing the radius of gyration or by in- 
creasing the wing to fuselage mass ratio. Damping of the torsion model is best 
when the center of gravity and elastic center are coincident, but damping of the 
bending mode is improved by moving the center of gravity and elastic center 
apart - also by decreasing the wing to fuselage mass ratio. 

As flight velocity increases the dynamic overstresses increase and the two 
oscillatory modes approach each other in frequency; each mode becoming more 
highly damped until one of the modes very suddenly becomes unstable and flutter 
occurs. Neglecting the unsteady forces yields a highly conservative estimate of 
flutter speed in the present case. 

Increasing altitude at constant velocity decreases the damping in each mode. 

This study has been concerned with changes in only one parameter at a time. 
Further stuffy might well be done on the interactions of two or more parameters: 
for example, the influence of the various parameters on flutter speed. Another 
problem of interest would be the behavior of the model under periodic gust loads. 

This procedure could be adapted to study of three-dimensional models by 
dividing the wing into segments. The feasibility of such a solution would be 
increased by use of the one-amplifier circuits for the unsteady forces. 
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Parameter 


Standard 

Value 


Range 


elastic center location 
(fraction of chord) 


.35 


. 25-. 50 


center of gravity location 
(fraction of chord) 


.45 


.30-. 50 


radius of gyration 
(fraction of chord) 


.25 


.20-. 40 


mass ratio, wing to 
fuselage 


.333 


. 15-. 75 


standard altitude 
(feet) 


S. L. 


0-30, 000 


flight velocity 
(ft/sec. ) 


380 


100-1200 


relative density ( y ) 


50 




wing chord (C) 

> 


11.9 




bending stiffness 
parameter (Kg) 


. 1 




torsion stiffness 
parameter (K/j.) 


.05 





Table I. Model Parameters 
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Distance Iron leading edge (chord lengths) 



Fig. 7. Change of Dynamic Overstresses in Response 
to Movement of the Center of Gravity or Elastic Center 
from standard position 
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Distance from leading edge (chord lengths) 

Fig. 8. Effect of Location of Center of 
Gravity and Elastic Center on Frequency 
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Fig. 9. Effect of Location of Center of 
Gravity and Elastic Center on Damping 
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Distance from leading edge (chord lengths) 



£Lg. 10. Effect of Location of Center of 
Gravity and Elastic Center on Mode Shape 
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Radius of gyration about the elastic center 
(chord lengths) 



Fig. ij. Effect of the Radius of 
Gyration on Dynamic Overstresses 
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Radius of gyration about the elastic center 
(chord lengths) 



Fig. 12 . Response of Frequency, Damping and 

Kode Shape to Change in Radius of Qyration 
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Fig. 1 3 . Effect of Hass Ratio on Dynamic Overstresses 
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Fig. 1 5. Effect of Hass Ratio on Mode Shapes 
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Fig. 1 6. Effect of Airspeed on .Dynamic Overstresses 
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Flight velocity (ft/sec) 



Fig. 17. Effect of Airspeed on Frequency and Damping 
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Fig. 18 . Effect of Airspeed on Mode Shape 
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Fig. 21 . Effect of Gust of Finite Gradient 
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